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Abstract 

Liquid foams have been observed to behave like immersed granular materials in at least one re- 
spect: deformation tends to raise their liquid contents, a phenomenon called dilatancy. We present 
a geometrical interpretation thereof in foams squeezed between two solid plates (2D GG foams), 
which contain pseudo Plateau borders along the plates, and in 3D foams. While experimental ob- 
servations evidenced the effect of a continuous deformation rate (dynamic dilatancy), the present 
argument applies primarily to elastic deformation (static dilatancy). We show that the negative 
dilatancy predicted by Weaire and Hutzler (Phil. Mag. 83 (2003) 2747) at very low liquid fractions 
is specific to ideal 2D foams and should not be observed in the dry limit of real 2D foams. 

PACS numbers: 83.80.Iz Emulsions and foams - 83. 60. He Normal stress differenees and their effeets - 
45.90.+t Other topics in classical mechanics of discrete systems 
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I. DILATANCY IN FOAMS AND GRAINS 

Liquid foams [ij] and granular materials are common examples of materials whose mechan- 
ical behaviours exhibit complex features, such as non-homogeneous flows: shear-banding, 
fracture or jamming. In the present work, we are interested in yet another such complex 
behaviour, called "dilatancy". Dilatancy was described by Bagnolds 2| in the context of 
granular materials: upon deformation, because grains are forced to move while avoiding 
each other, the medium swells to some extent. In other words, the fluid volume fraction 
is increased. This effect can remain unnoticed in air. By contrast, a spectacular absorption 
of liquid is obtained upon deformation of an immersed granular sample. 

In liquid foams, it is not a priori so obvious whether dilatancy should be expected or not. 
Because bubbles are deformable individually, a foam should be able to deform substantially 
without altering its liquid fraction. 

In 3D foams, dynamic dilatancy was indeed observed under stationary shear j^: in a 
device where the shear rate localizes, the liquid volume fraction was observed, both visually 
and through electrical conductimetry, to stabilize, within seconds, at a higher value in 
the region being sheared continuously, while it remained lower in both statically deformed 
regions (and when shearing was stopped, the liquid fraction was observed to become uniform 
again within seconds). Dynamic dilatancy, whose microscopic origin is yet uncertain 18| . 
thus constitutes one out of many mechanisms that could favour shear banding, a family of 
rheological behaviours commonly observed in many other complex fluids such as giant micelle 
solutions, multilamellar vesicles or onions, granular pastes, etc. The stabilizing mechanism 
is as follows. As both the foam shear modulus and its yield strain are known to decrease 
as the liquid fraction increases a region that contains more liquid will host statistically 
more plastic events than other regions and, if dynamic dilatancy is indeed present, even 
more liquid will then permeate towards that region from less active regions. This process 
reaches a stationary regime (co-existing shear bands) when the dilatancy-induced pressure 
gradient is balanced by the osmotic gradient associated with fluid concentration gradients. 

Whether static dilatancy (resulting from a fixed deformation) also exists in liquid foams 
has not been determined experimentally so far. It might constitute a destabilizing factor for 
a homogeneous flow and generate shear banding in the first place. 

A theory of static dilatancy was derived from thermodynamical considerations a few 
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FIG. 1: Four types of 2D foams, (a) Top- view of an ideal, truly two-dimensional foam: either 
the third dimension is purely absent, or the foam is invariant with respect to the third dimension, 
(b) Side-view of a monolayer of bubbles between two glass plates ("GG" foam), (c) Monolayer 
of bubbles floating on a liquid bath below a glass plate ("LG" foam), (d) Monolayer of bubbles 
floating on the free surface of the liquidhath, in contact with air ("LA" foam). Dilatancy is no real 
issue for LG or LA foams because of the underlying reservoir of liquid. By contrast, the GG foam 
holds liquid not only in the Plateau borders just like an ideal 2D foam (second term of Eq. ([I|)), 
but also in the pseudo Plateau borders which run along the edge between two bubbles and touch 
one solid plate each (first term of Eq. 

years ago P, 0] for an ideal, truly two-dimensional foam (Fig. [1^), with predictions based 
on empirical or numerical results concerning the osmotic pressure and the elastic modulus. 
The model leads to the surprising prediction that the sign of the dilatancy effect should be 
negative for very low liquid fractions, as a consequence of the fact that the elastic modulus 
does not depend on the liquid fraction in the dry limit (Decoration Theorem l|, |6|). 

In this paper, we provide a simple geometrical description of static dilatancy in a dry 
liquid foam. An elastic deformation of such a material induces an elongation of the bubbles 
while conserving their volume constant, thereby increasing the total Plateau border lengh 
per bubble, or per unit volume. Since most of the liquid is located in Plateau borders, the 
liquid fraction is thus increased. 

According to this description, we calculate the magnitude of the dilatancy effect analyt- 
ically for a crystalline, hexagonal GG foam (Fig. [lb) as well as for a 3D Kelvin foam, and 
we discuss negative dilatancy 
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II. STATIC DILATANCY IN A DRY 2D GLASS-GLASS FOAM 



Here, we introduce a very simple geometrical description of static dilatancy. In a "GG" 
foam, there are two contributions to the liquid content, which are presented on Fig. [D 

^^Hq ^ P R% + (2v^ - tt) R^H (1) 

where finq is the volume of liquid per bubble. The first term corresponds to the pseudo 
Plateau borders (with the corresponding radius of curvature -Rps and the bubble perimeter 
P) while the second term corresponds to the genuine Plateau borders (radius R and gap H 
between the solid plates). 

In the present work, we will assume both P ^ R {i.e., (j) <^ 1) and -Rps = R H (floor 
tile regime). In other words, we concentrate on regimes E, F and G of Fig. HI Regime E 
corresponds to a floor tile GG-foam, while regime G is the truly 2D limit of refs. 0, a| 

We here focus on floor tile "GG" foams (regime E of Fig. H]), i.e., with not only a low 
liquid volume fraction (0^1, achieved when R <^ P) but also small pseudo Plateau 
borders compared to the gap between the solid plates (-Rps ^ H), in which case one has 
-Rps ^ R. The opposite limit of "pancakes" (-Rps ^ H/2 < R, regime ABCD of Fig. |4j) will 
be discussed elsewhere 

Let us consider such a 2D glass-glass foam free of any significant anisotropic {i.e., devia- 
toric) in-plane stress. If the foam is now deformed elastically, the average bubble perimeter 
P increases as the typical bubble elongates while conserving its volume constant. Since the 
total amount of liquid per bubble remains constant on short time scales, the increase in 
perimeter causes the pseudo Plateau borders to shrink accordingly (see Fig. [2]). Later on, 
as permeation takes place, the pseudo Plateau border radius -Rps may return to its original 
value, and the part of the liquid that is located in the pseudo Plateau borders increases 
accordingly. 

Thus, soon after deformation, the pressure difference between the gas and the liquid 
must increase as the pseudo Plateau borders (and Plateau borders) shrink. Later, after 
permeation has taken place, the liquid content of the foam has increased as compared to its 
value when the foam was at rest. In other words, there is a positive dilatancy effect as soon 
as the perimeter increases. 

In order to evaluate the change in bubble perimeter, let us consider a dry, crystalline 
foam subjected to an arbitrary elastic, homogeneous deformation. Up to a global rotation, 
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FIG. 2: Pseudo Plateau borders between two bubbles squeezed between two glass plates. Original 
configuration, with top- view (al) and perspective (a2). When a bubble which is initially at rest is 
suddenly deformed, its perimeter, i.e., the total length of its pseudo Plateau borders, is increased. 
On short time scales, the amount of liquid that surrounds the bubble remains constant, hence 
the quantity of liquid per unit length of the perimeter decreases (bl). Thus, the pseudo Plateau 
borders shrink, and the corresponding radius of curvature Rps is reduced (b2). This implies that 
the pressure in the liquid is lowered. As a result, liquid from less deformed regions of the foam, 
where the liquid pressure is not as low, can soon permeate and compensate for the Plateau border 
shrinking (c). 

this deformation is equivalent to an elongation by some factor A in one direction and to a 
compression by the same factor in the perpendicular direction. In such a hexagonal foam, 
in the dry limit, the bubble perimeter increases in the following way: 
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The perimeter thus increases when the foam is deformed (A 7^ 1), as can be seen on Fig. [31 
This law, which is exact in the dry limit, holds as long as no Tl process occurs and does 
not depend on the orientation of the crystalline network with respect to the direction of 
elongation. In fact, it even holds for polydisperse hexagonal foams j^. The elongation at 
which Tl processes occur depends both on orientation or on liquid volume fraction. In the 
dry limit, it is A = y/S when elongation is perpendicular to some facets and larger for other 
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elongations. This yield deformation is lowered as the foam becomes wetter as described 
already long ago The expected perimeter increase for A = is of order 15%, which 
means that gas-liquid pressure difference is expected to increase by 7% on short time scales, 
while the foam liquid fraction is expected to increase by 15% after permeation. 



III. DILATANCY IN THREE DIMENSIONS 

In a dry, three-dimensional foam, most of the liquid is located in the Plateau borders. 
Following the same line of thought as above, static dilatancy is expected if deformation 
induces an increase in the total Plateau border length per bubble (or per unit volume of 
foam) . 

To test this in a simple manner, we took an approximate version of Kelvin's bcc bubble 
packing. Kelvin's packing is one of the most commonly encountered in monodisperse 3D 
foams, even though it does not have the lowest energy . In a real Kelvin foam, square faces 
have edges bending slightly outwards, while hexagonal faces are slightly non-planar. Here, 
each bubble is assumed to be a truncated octahedron, called tetrakaidecahedron, with six 
planar, square faces and eight regular hexagonal, planar faces. This approximation was used 
and discussed by Reinelt and Kraynik 



lo| to express the stress and the energy under large 
elastic deformations. Here, we present the changes in the total Plateau border perimeter. 

Computing the bubble perimeter (sum of all edge lengths) as a function of elongation 
has been performed for three orientations of the elongation (with a compression by I/VA 
in the perpendicular directions. As can be seen on Fig. [3l the result is very similar to the 
two-dimensional case, and indicates that significant dilatancy is to be expected in dry 3D 
foams, too. Incidentally, Fig. [3] indicates that for the present approximation of a Kelvin 
foam, the response is independent of the direction of elongation: 

Pf{X) _ P|°(A) _ P^J^iX) _ 1 , 2 ] 



P3D p3D p3D 3 



IV. OSMOTIC PRESSURE AND DILATANCY 



A+ ^ 



The dilatancy coefficient is defined as the (second) derivative of the osmotic pressure 
with respect to the deformation e: x = d'^'^osm/de'^ . With the variations of the average 
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FIG. 3: Sum of all Plateau border lengths in a 2D (hexagonal, Eq. dJ])) foam or 3D (Kelvin, Eq. ([3])) 
foam, as a function of the elongation A. In 3D case, elongation is performed in direction A (normal 
to a square facette), B (normal to a hexagonal facette) or C (one of the edges of Kelvin's cell) with 
compression l/^fX in both perpendicular directions. 

bubble perimeter in mind, let us now derive the variation of the osmotic pressure in a floor 
tile GG-foam. From Eq. ([2]), expressed in terms of A = 1 + e, one has bPjP ^ Hence, 
X = SiTosm/iSP/P), where the variation ^VTosm is taken at constant liquid volume fraction 0. 

When the foam is conflned in a container, the osmotic pressure TTosm corresponds to the 
force that must be applied externally to one of the conflning walls if the latter is permeable 



to the liquid but not to the bubbles. ' 

r 

tensor and not just a scalar quantity [4 



le osmotic pressure (which is in fact a symmetric 
^) is thus the difference between the stress in the 



foam and the pressure applied by the pure liquid on the other side of the semi-permeable 
wall: TTosm = —pi — 0" (where tensile stresses and pressures are both counted positively). The 
stress in the foam includes a pressure contribution from the liquid (pi) and from the gas (pg), 
as well as a tensile contribution from the interfaces: a = —(ppil — (1 — <P)PgI + o"™*'^''^. In the 
floor tile limit {R <^ H), the interface stress, averaged over in-plane orientations and over the 
sample thickness, is related to the bubble perimeter [3]: < cr™*°''^ >2D— (1 ~ 0) ["^ifr] + 
where 7 is the surface tension and Vt is the bubble volume. Hence, osmotic pressure (now 
expressed as a scalar quantity) can be written as: < vTosm >2D= (1 — 4>) [Ap — -^^f-] — 
where Ap = Pg — p\ = = ^ is the pressure difference between the gas and the liquid. 
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FIG. 4: Regimes of a GG-foam with low liquid fraction {(f> <^ 1), in terms of the bubble perimeter 



P (measured on top view), the Plateau border radius R (as seen from above too) and the cell heig^ 



ht 



H. Real GG-foams {P > H) can be found in two main configurations: the "pancake" regime 
{H - 2Rps <^ H <^ R, region ABCD) and the "floor tile" regime E {R = Rps < H). Regimes 
F and G do not correspond to real foams (which cannot sustain the condition H ^ P without 
destabilizing into 3D foams), but regime G is the limit of truly two-dimensional foams (H oo) 



addressed in earlier works on dilatancy 
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The dilatancy then results from Eqs. ([T]) and 

S < VTosm >2D 

X = 



SP/P 



i/m jp 



H 2A 

l-7r/4 P 



The last term, in which A= Q/H, originates in the interface stress contribution (j™*^^ and 
is equal to one half of the interfacial energy per unit height (7P per bubble). 



V. NEGATIVE DILATANCY? 



For a GG-foam (regime E with P ^ H, see Fig. H]), one has x — Aj9/2 > 0. By contrast, 
when H ^ P (regime G), one gets the limit of a truly two-dimensional foam since the 
pseudo Plateau borders (total length P per bubble) are too small compared to the Plateau 
borders (length H) to contribute to dilatancy significantly. Hence only the last term remains 
in Eq dH), and dilatancy is negative: x — ~'2j- ^^^^ recover exactly the result x — — 2G 

9 



by Weaire and Hutzler {4] in the dry limit [19|]. For instance, in the case of a crystalhne, 
hexagonal 2D foam: x — — 4v^7/-P- 



VI. DISCUSSION 



In the present work, we have shown that positive static dilatancy in liquid foams might 
result from the overall increase in Plateau border lengths when the foam is being deformed 
elastically. In 2D GG foams, this would be in fact the pseudo Plateau borders (those 
along the solid walls). The present study was restricted to crystalline foams, whether 2D 
GG hexagonal foams or 3D Kelvin foams. A weaker (or even reversed) effect cannot be 
excluded, at this stage, in the case of disordered or polydisperse foams, or even simply for 
different 3D crystalline structures. Indeed, in 3D at least, there seems to be no trustworthy 
reason why the minimum in the total facette surface area (which is the definition of a foam 
at rest, in the dry limit) should coincide with the minimum in the sum of all Plateau border 
lengths. For a GG foam in a somewhat less dry configuration (regimes A-D of Fig. HI), the 
expected dilatancy effect can be shown to be less intense l7|. 

The negative dilatancy predicted [4, |5| in the dry limit in truly 2D foams (invariant along 
the third dimension of space) coincides with one particular limit of our study, which indeed 
corresponds to a truly 2D foam. It is then related to the interfacial (and deviatoric) part 
of the stress. By contrast, in the other regimes, we obtain a positive dilatancy, related to 
the pressure in both fluid phases. Let us mention that negative dilatancy is also possibly 
expected in the wet limit of the pancake regime [3]. 

If static dilatancy is conflrmed experimentally, it should be coupled with drainage models 



ahich describe permeation) and incorporated into rheological models of foams 
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Eq. ^ yields 6P/P = F^/S. Hence, G = \ -fP/A and x 
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